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Networks are topological and geometric structures used to describe systems as different as the In¬ 
ternet, the brain or the quantum structure of space-time. Here we define complex quantum network 
geometries, describing the underlying structure of growing simplicial 2-complexes, i.e. simplicial 
complexes formed by triangles. These networks are geometric networks with energies of the links 
that grow according to a non-equilibrium dynamics. The evolution in time of the geometric net¬ 
works is a classical evolution describing a given path of a path integral defining the evolution of 
quantum network states. The quantum network states are characterized by quantum occupation 
numbers that can be mapped respectively to the nodes, links, and triangles incident to each link of 
the network. We call the geometric networks describing the evolution of quantum network states the 
quantum geometric networks. The quantum geometric networks have many properties common to 
complex networks including small-world property, high clustering coefficient, high modularity, scale- 
free degree distribution. Moreover they can be distinguished between the Fermi-Dirac Network and 
the Bose-Einstein Network obeying respectively the Fermi-Dirac and Bose-Einstein statistics. We 
show that these networks can undergo structural phase transitions where the geometrical properties 
of the networks change drastically. Finally we comment on the relation between Quantum Complex 
Network Geometries, spin networks and triangulations. 
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I. INTRODUCTION 

Networks are discrete structures that can be used to 
describe, model and understand a variety of real systems, 
including complex interacting systems w or the micro¬ 
scopic nature of space-time [sm. 

Recently, in network science the characterization of the 
complexity of networks with geometrical and topological 
methods is gaining large momentum with several works 
related to the definition of curvature of the networks [IQ^ 
[22], persistent homology [23H^ . complex networks em¬ 
bedded in finite dimensions [26U29] , and the study of the 
hyperbolicity of complex networks j3QH38] . 

In this context, it is becoming clear that in order to 
characterize the geometry of networks it is important to 
describe the underlying structure of simplicial complexes. 
A simplicial complex is constructed by gluing together 
simplices such as points, lines, triangles etc., along their 
faces. Therefore different works characterize ensembles of 
random simplicial complexes [39ti4^ and their geomet¬ 
rical and topological properties. Recently a model for 
emergent complex network geometry has been proposed 
based on growing simplicial complexes [43] . 

In quantum gravity a central problem is to find the ap¬ 
propriate model describing the geometry of space-time 
at the quantum level. Different approaches have been 
proposed in which geometry emerges from some prege¬ 
ometric phase [521154] that include spin networks and 
loop quantum gravity ISE], spin foams [44|, causal dy¬ 


namical triangulations gaisn], causal sets m, energetic 
causal sets [46H48] , network cosmology m, and quantum 
graphity ISSHSU- Network-like structures play a funda¬ 
mental role in all these approaches. 

This suggests that the emergence of geometric struc¬ 
ture from the quantum description of networks is a more 
general mathematical problem that can be not only rel¬ 
evant for understanding the structure of space-time, but 
which might also help to understand general complex net¬ 
work structures. 

Already in the early days of the field of network sci¬ 
ence the relation between complex network topologies 
and quantum statistics was shown in the framework of 
the Bianconi-Barabasi model (58] [59| that is a growing 
network model with preferential attachment and fitness 
of the nodes which display a Bose-Einstein condensation. 
This model can also be extended to weighted networks 
[62] described by the Bose-Einstein statistics and under¬ 
going also the condensation of the weight of the links. 
The relation between growing Cayley trees with fitness 
of the nodes and Fermi-Dirac statistics has been found 
in m and the underlying symmetries between the mod¬ 
els in m and m have been discussed in m- In the 
context of equilibrium network models it has been shown 
that quantum statistics emerges to describe simple or 
weighted networks [63] . 

Here we characterize the non-equilibrium evolution of 
networks constructed from growing simplicial complexes 
of dimension two, i.e. formed by triangles, and such that 


2 


to each link we associate an energy e. We show that ge¬ 
ometrical complex networks emerge from this dynamical 
evolution which display at the same time small-world net¬ 
work properties [70] , exponential or scale-free degree dis¬ 
tribution EB, high clustering coefficient and high mod¬ 
ularity [5]. These networks can be either planar or non- 
planar with an Euler characteristic that is either x = 1 
(planar) or x oc A/", where N is the network size, indi¬ 
cating a finite average curvature in the network. As we 
will show in two limiting cases of this network dynamics, 
these networks describe the evolution of quantum net¬ 
work states. These network states are constructed along 
similar lines used in the quantum gravity literature [SSl- 
EB by associating an Hilbert space to each node of the 
network and two Hilbert spaces to each possible link of 
the network. These network states evolve through a non 
equilibrium, Markovian dynamics. The network states 
can be mapped to geometric networks and have an evo¬ 
lution described by an appropriate path integral. The 
network dynamics describes the paths of single histories 
of networks on which the path integral is calculated. 

We distinguish between Fermi-Dirac Networks and 
Bose-Einstein Networks. For both of them the number 
of triangles incident to a given link is n + 1. However, for 
Fermi-Dirac Networks n can only take the values n = 0,1, 
whereas for Bose-Einstein Networks n can take any in¬ 
teger value n = 0,1, 2 ... These networks evolve in such 
a way that at the global scale the average of n over the 
links with energy e follows the Fermi-Dirac statistics for 
the Fermi-Dirac Network and the Bose-Einstein statistics 
for the Bose-Einstein Network. 

These network structures depend on an external pa¬ 
rameter that we call the inverse temperature p. As a 
function of P they undergo major structural phase tran¬ 
sitions in which the network structure changes drasti¬ 
cally. In the case of the Fermi-Dirac Network, for f3 > f3c 
the network is not any more small world, but acquires a 
finite Hausdorff dimensionality. In the case of the Bose- 
Einstein Network, for > /3c a link acquires a finite 
fraction of triangles, and the nodes at the end of the link 
acquire a finite fraction of all the links. This geometrical 
phenomenon is the Bose-Einstein condensation for these 
networks. 

We observe here that the quantum network states stud¬ 
ied in this paper are by no means the only way to asso¬ 
ciate a quantum state to a network. In particular in 
the quantum computation community, alternative ap¬ 
proaches [64II68] have been widely explored, characteriz¬ 
ing quantum transport, quantum random networks, and 
quantum networks in which the links correspond to en¬ 
tangled states. 

The paper is organized as follows. In Sec. H we de¬ 
scribe the geometric network model with energy of the 
links depending on the parameter m fixing the maximum 
number of triangles incident to a link. Moreover we de¬ 
fine the entropy rate of the model and we describe the 
observed phase transition. In Sec. HI we define the quan¬ 
tum network states. In Sec. IV. we define the evolution 


of the Fermi-Dirac quantum network state and the Bose- 
Einstein quantum network state. In Sec. V we study the 
Fermi-Dirac Network (given by the geometric network 
model with m = 2). We show that this network charac¬ 
terizes the Fermi-Dirac quantum state, and we show that 
it is globally described by the Fermi-Dirac statistics. Fi¬ 
nally we compare the analytical results to simulations 
and we describe the phase transition occurring at low 
temperatures. In Sec. VI we study the Bose-Einstein 
Network (given by the geometric network model with 
m = oc) showing that it fully characterizes the Bose- 
Einstein quantum state, follows the Bose-Einstein statis¬ 
tics and undergoes the Bose-Einstein condensation at low 
temperatures. In Sec. VH we consider the thermody¬ 
namics of the networks and we consider the case in which 
we project the quantum network state on an unlabeled 
final network state. In Sec. VHI we generalize the geo¬ 
metric model introducing a new parameter p and a new 
process of addition of triangles that allows for the gen¬ 
eration of network geometries that are not planar. We 
characterize the geometry of these networks and we de¬ 
scribe the phase transitions observed at low temperature. 
In Sec. IX we generalize the evolution of the quantum 
network states corresponding to the generalized geomet¬ 
ric network model. In Sec. X we describe the dual of the 
networks generated by the proposed model and we com¬ 
ment on the relation between the Fermi-Dirac Network 
and spin networks. In Sec. XI we comment on the re¬ 
lation between Complex Quantum Network Geometries, 
triangulations and foams. Finally in Sec. XH we give the 
conclusions. 


II. GEOMETRIC NETWORK WITH ENERGY 
OF THE LINKS 

A. Evolution of the geometric networks with 
energy of the links 

Real networks display at the same time several struc¬ 
tural properties (including finite clustering coefficient, 
significant modularity, finite spectral dimension, hetero¬ 
geneous degree distribution) that have been shown to 
be captured by a very simple model of emergent geome¬ 
try recently introduced by the authors [43]. The model 
proposed in m is a non-equilibrium model of growing 
simplicial complexes of dimension = 2, i.e. formed 
by gluing triangles along their edges. In this model each 
link can belong at most to a number m of triangles where 
the parameter m can take any finite value m > 2 or the 
value m = oc, indicating the case in which each link can 
belong to an arbitrarily large number of triangles. In 
the case m = 2 the model reproduces random manifolds 
of dimension = 2 with an exponential degree distri¬ 
bution and random distribution of local curvatures, in 
the case m = oc the model generates scale-free networks 
with finite clustering coefficient and significant modular¬ 
ity quantifying the relevance of their community struc- 
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ture. 

In [43] all the nodes and all the links are treated 
equally, having the same probability to attract new tri¬ 
angles. Nevertheless, in complex systems, attaching a 
new triangle to a given link might not have the same 
probability of attaching it to another link. 

Already in the context of complex networks growing 
by preferential attachment w, the heterogeneity of the 
nodes in attracting new links has been recognized to be 
essential to characterize the evolution of networks, as for 
example the World-Wide-Web or the Internet [58l l59] . 
Usually, this heterogeneous ’’quality” of the nodes is 
modeled by associating each node to an energy drawn 
from a given distribution. Interestingly, complex net¬ 
works with energy of the nodes have been shown [59l |60| 
to be characterized by quantum Bose-Einstein and Fermi- 
Dirac statistics, and might display a Bose-Einstein con¬ 
densation in which one node grabs a finite fraction of 
the links. This phase transition is relevant for a number 
of complex networks including economical, technological 
and social networks in which nodes connected to a finite 
fraction of the nodes might emerge. 

In the quantum gravity literature, the relation between 
networks and quantum states has been recently explored 
[551157] to construct models of emergent space-time ge¬ 
ometry. In these works, each network is associated to 
a quantum network state and the network structure is 
dictated by an equilibrium Hamiltonian dynamics. 

Here we consider networks constructed by a non- 
equilibrium dynamics describing the underlying struc¬ 
ture of simplicial complexes constructed by the addition 
of connected complexes of dimension = 2, i.e. tri¬ 
angles. These networks display non trivial geometrical 
properties, characterizing in some limit planar random 
manifolds, as will be discussed later in the paper. For 
this reason we call them geometric networks. As in the 
geometric network model [43] we assume that each link 
can belong at most to a number m of triangles. Moreover 
we associate energies both to nodes and links describing 
the different ability of nodes and links to attract new 
triangles. 

In studying this model our goal is two-fold. On one side 
we aim at characterizing a wider class of emergent ge¬ 
ometries, and their possible structural phase transitions 
in order to unveil the basic geometric properties of com¬ 
plex networks. On the other side we aim at furthering 
our understanding on the relation between the network 
evolution and quantum mechanics by exploring the con¬ 
nection between network evolution, quantum statistics, 
and evolution of quantum network states constructed us¬ 
ing methods similar to the one introduced in [551157] . 

The energies of the nodes and of the links are defined 
as follows. Every node i of the network is associated with 
the energy of the link uji > 0 drawn from a distribution 
g{uj) [69]. The energy uJi is assigned to the node i when 
the node is added to the network and is quenched during 
the growth of the network. Every link £ = (i, j) between 
node i and node j is associated with the energy of the link 


Cij which is a given symmetric function of the energy of 
the two nodes i and j, i.e. 

<^ij = ( 1 ) 

with eij > 0. 

We define the so called spin Jij of the link = (i, j) as 

~ 2 ^ ) • (^) 

The spins of the links belonging to a triangle between the 
nodes i, j and r satisfy the conditions 

\Jir '7jfr| ^ Jij ^ Jir ~£' Jjr- (3) 

This result remains valid for any permutation of the or¬ 
der of the nodes and r belonging to the triangle. Al¬ 
though most of the derivations shown in this paper can 
be performed similarly for either continuous or discrete 
energy of the nodes and of the links, here we consider 
the case in which the energies of the nodes {uJi} and the 
energy of the links {e^j} are discrete. In particular, if 
the energy of the nodes takes integer values, the spin 
of the links takes half-integer values and Eqs. Q can 
be interpreted as the Clebsch-Gordon relations between 
the half-integer spins of the links of each triangle. This 
property motivates dubbing this variable a spin. 

Specific expressions of the energy eij of the link {i^j) 
might depend on the spin Jij of the link. Examples of 
specific choices for the energy of the link are the quadratic 
relation, 

^ij = Jij {Jij + 1) (4) 

or the linear relation 

^ij — ‘^Jij — ^i T • (5) 

Here we want to keep the generality of the model and we 
will take eij given by Eq. Q unless a specific functional 
form of the energy of the link is indicated. 

The geometric network model is the underlying net¬ 
work of a simplicial complex of dimension ^2 = 2 formed 
by gluing triangles along the edges. We assume that each 
link can belong at most to a number m of triangles where 
the parameter m can take any finite value m > 2 or the 
value m = oc, indicating the case in which each link can 
belong to an arbitrarily large number of triangles. We 
call the links to which we can still add at least one trian¬ 
gle unsaturated. All the other links we call saturated. In 
the case m = oc, all the links are unsaturated. We start 
at time t = 1 from a network formed by a single triangle, 
a simplex of dimension dn = 2. At each time we add a 
triangle to an unsaturated link (i^j) of the network. We 
choose this link with probability given by 

tt[i] _ + 

~ Z 


( 6 ) 
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where Z = Z is given by 

Z = ^ ^ e ^ + Tir/^g^. (7) 

r<s 

Here we introduced several time-dependent quantities 
which we will use repeatedly in this paper: a^j is the 
element (i^j) of the adjacency matrix a of the network, 
^ij is equal to one (i.e. ^ij = 1) if the number of triangles 
to which the link (i, j) belongs is less than m, otherwise it 
is zero (i.e. ^ij = 0), riij + 1 is equal to the total number 
of triangles incident to the link (i, j). Having chosen the 
link (i,j) the simplicial complex at time t is constructed 
by adding a node r, two links (i,r) and (j, r) and the 
new triangle linking node i, node j and node r. The ge¬ 
ometric complex network is the network structure of the 
resulting simplicial complex. 

Therefore the number of nodes N of the network grows 
linearly with time and is given hy N = t 2. 

The linking probability depends on the parameter (3 > 
0 that we call inverse temperature. For = 0, all the 
links that are unsaturated have equal probability to be 
selected. For /3 > 0 instead, unsaturated links with low 
energy eij are more likely to be selected than links with 
higher energy. 

With the above algorithm we describe a growing sim¬ 
plicial complex formed by adding triangles. From this 
structure we can extract the corresponding network 
where we consider only the information about node con¬ 
nectivity (which node is linked to which other node). We 
call this network model the geometrical growing network. 
In Figure we show schematically the dynamical rules 
for building the growing simplicial complexes and the 
growing geometrical networks that describe its underly¬ 
ing network structure. In the following we will focus 
in particular on the limiting cases in which m = 2 (2- 
dimensional manifolds), or m = oo. For reasons that will 
become clear in the following, we will call the growing 
geometric network with m = 2 the Fermi-Dirae Network 
and the one with m = oc the Bose-Einstein Network. 
In Figure we show examples of the first few steps of 
their evolution. The Fermi-Dirae Network and the Bose- 
Einstein Network will also be indicated as quantum geo¬ 
metric networks. 


B. Entropy rate of the network evolution 

Entropy measures for network evolution are very im¬ 
portant characteristics for evaluating the interplay be¬ 
tween randomness and order in these structures [zacs]. 
In particular, for growing network models, the entropy 
rate m characterizes how the space of typical network 
dynamical evolutions increases with time. A change in 
the scaling of the entropy rate typically indicates a phase 
transition in the network m- The geometric network 
evolution is described by the sequence {uj{t')^l{t')}t'<t 
where uj(t) indicates the energy of the node added to the 
network at time t, l{t) = (i^j) indicates the link chosen 
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FIG. 1: (Color online)The growing geometrical network is the 
underlying network structure of a growing simplicial complex 
in which triangles are continuously attached to the simplicial 
complex and glued to one unsaturated link. The link where 
the new triangle is added is chosen with probability 
given by Eq. (^. The figure shows an example where the 
maximum number of triangles incident to a link is m = 2. 


at time t with probability given by Eq. (j^. At any 

given time, therefore, £{t) indicates the link to which the 
new triangle is attached. 

The entropy rate of the network evolution can be ex¬ 
pressed as 

Hoit) = - 

X (8) 

where P{uj{t),i{t)\{uj{t')^i{t')}t'<:t) is the probability 
that, given the temporal evolution of the network un¬ 
til time t — 1, at time t a new triangle is attached to the 
link £{t) with the new node of this triangle having energy 

At time t the probability that the new node has en¬ 
ergy uj{t) = uj is given by the probability distribution 
^(cj), and it is independent of the previous evolution of 
the network. The probability of choosing the link £{t) is 
given by that depends on the previous history of 

the network. Moreover uj{t) and £{t) are independent. 
Therefore the entropy rate can be written as 

Hoit) = ( 9 ) 

with specified below. In particular is 

the contribution to the entropy rate due to the random 
distribution of the energy of the nodes, it is independent 
of time and is given by 

=-y]s((w)ln5(w). (10) 

UJ 

The quantity of the growing geometric network 

evolution defines the contribution to the entropy rate 
due to the choice of the link where the new triangle is 
attached and is given by 

«'■'(') = <“) 

l<j 

evaluated at time t. With Eq. ([^, we get 
ffW(i) = /3{eij)+\ogZ 


( 12 ) 
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FIG. 2: (Color online) The Fermi-Dirac Network evolution 
(panel A) and the Bose-Einstein Network evolution (panel B). 
At each time a new triangle is added to a link (z,B chosen 
according to the probability given by Eq. JgI. The 

maximum number of triangles incident to a link is m = 2 
for the Eermi-Dirac Network evolution and m — oo for the 
Bose-Einstein Network evolution. 


where 

i<j 


(13) 


We note here that, as the inverse temperature P changes, 
we might expect a phase transition in the network char¬ 
acterized by a different scaling of the entropy rate 
and the normalization constant Z with time t below and 
above the transition. 

The normalization constant Z is fixed by Eq. 0. For 

= 0, Z grows linearly with time t. In fact for /3 = 0 and 
m = 2, riij = 0 only if ^ij = 1 and riij = 1 only if ^ij = 0. 
Moreover, since at each time we add two unsaturated 
links and we remove one unsaturated link. 


Z — + riij) — t + 2. (14) 

i<j 

For P = 0 and m = oo instead, all the links are unsatu¬ 
rated, i.e. ^ij = 1 and every triangle is incident to three 


links. Therefore, since we add a triangle for every time 
step. 


Z = + Uij) = 3t. (15) 

i<j 

For a significant range of values of /3 > 0 we will still 
continue to have Z (x t for t ^ 1 because Z is a sum over 
a linearly growing set of non-zero variables. For /3 ^ oo 
however, only the unsaturated links with minimal energy 
of the links €ij = cq will contribute to the sum defined in 
Eq. 0 because the dynamics becomes extremal. There¬ 
fore we can have Z 0{1). In this case a phase transi¬ 
tion is expected to occur in the network at the structural 
level. This phase transition induces a substantial change 
in the geometry of the networks above and below the 
transition as discussed in the next paragraph. 


C. Phase transition in geometric networks 

The networks constructed according to the model de¬ 
fined in Sec. |II A| are planar. In fact the Euler number 
of the simplicial complex from which they are extracted 
is constant during the network evolution and is given by 

X = N - L^T = 1 (16) 

where N is the total number of nodes, L is the total 
number of links and T is the total number of triangles. 

Let us prove this result recursively. At time t = 1 
the geometrical networks are formed by a single triangle 
A/' = 3,I/ = 3,T = 1. Therefore we have 

X = 1 - ( 17 ) 

At each time step we add a single node, two links and 
one triangle, therefore 

Ax = x{t) - x{t - 1) = 0. (18) 

This shows that the simplicial complexes constructed by 
gluing new triangles to a single existing link of the net¬ 
work have Euler characteristic % = 1. When consider¬ 
ing the network underlying each of these simplicial com¬ 
plexes, and embedding it in a plane, one can see that the 
number of faces F of the embedded graph is in fact equal 
to the number of triangles T of the simplicial complex 
when we do not count the external face of the planar 
network. Therefore these networks are planar. 

Another, equivalent way to prove that our networks are 
planar is to observe that these networks, by construction, 
do not contain any complete graphs of five nodes (sub¬ 
graph iFs) or any bipartite complete graph of six nodes 
(subgraph 

Additionally we define the boundary of the network, 
as the set of unsaturated nodes and links. Unsatu¬ 
rated links, are links (i, j) with ^ij = 0, while unsatu¬ 
rated nodes are nodes with at least an incident unsatu¬ 
rated link. Note that in the geometric networks studied 
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here, all the nodes are unsaturated since by construc¬ 
tion they are always incident to exactly two unsaturated 
links. Therefore all the nodes of the network belong to its 
boundary. For these networks the curvature Ri min] 
associated to each node i is given by 


ki Ti 4: — ki 
2 3 6 


3-T, 

6 


(19) 


where ki indicates the degree of node Ti indicates the 
total number of triangles incident to node i, and where 
the last equation can be derived by considering that in 
the present model Ti = ki — 1 for every i. The last expres¬ 
sion in Eq. (19), relating the curvature Ri of node i to 


the number Ti of triangles incident to it, has an intuitive 
explanation. As all triangles are isosceles, and each node 
is at the boundary of the network,each node incident to 
exactly = 3 triangles will have zero curvature, since 
the sum of the angles incident to it is tt. 

Here we focus on the quantum geometric networks 
(cases m = 2 and m = oc) and we study the geome¬ 
try of these network models as a function of (3. These 
networks are generated by a non-equilibrium dynamics 
that does not contain any indication about any embed¬ 
ding space. In the case m = 2, the Fermi-Dirac Neworks 
are planar manifolds describing random geometries. In 
the case m = oc, the Bose-Einstein Networks are planar 
scale-free networks but are not manifolds. 

Here we show numerical evidence that for given dis¬ 
tribution g{uj) and energy of the links eij = f{uji,ujj) a 
structural phase transition can occur in quantum geomet¬ 
ric networks. Specifically, we consider the case in which 
uj can only take integer values and the distribution g{uj) 
is Poisson with average c, i.e. 


UJ\ 


( 20 ) 


Moreover we take the energy eij of the generic link (i, j) 
given by Eq. (§. 

As a function of /3 we observe a phase transition in both 
the Fermi-Dirac Network and the Bose-Einstein Network. 
For [3 > Pc the structure of the network and its geometry 
change drastically as can already be seen from the visu¬ 
alizations of the networks (Figure]^ for the Fermi-Dirac 
Network and Figure for the Bose-Einstein Network). 
The transition is characterized by a different scaling of 
the entropy rate below and above pc- For P < Pc 
increases with time as IrL(t), due to the lin¬ 

ear scaling of Z oc t, while, for P > Pc^ = 0{1) and 
fluctuates widely during the network evolution. Here we 
discuss in detail the consequences of this transition in the 
Fermi-Dirac Network and in the Bose-Einstein Network. 
In Figure we show major geometrical and structural 
properties of the network as a function of the inverse 
temperature P across the phase transitions. In particu¬ 
lar we display the maximal shortest (hopping) distance 
from a given node of the initial triangle D, the maxi¬ 
mal degree kmax of the network, the entropy rate 


the modularity M US] calculated using the Louvain algo¬ 
rithm HU and the average clustering coefficient C across 
the phase transitions. In Figure we show major ge¬ 
ometrical and structural properties of the network as a 
function of time for given values of the inverse tempera¬ 
ture P below and above the phase transition. Finally in 
Figure and Figure we show the degree distribution 
P(/c), the average clustering coefficient C{k) of nodes of 
degree /c, and the distribution of the curvature P{R), for 
the Fermi-Dirac and the Bose-Einstein Network above 
and below the phase transition. 

For the Fermi-Dirac Network, the most important indi¬ 
cator of the phase transition is D which grows logarithmi¬ 
cally with time for p < pc^ and as a power-law for p > pc. 
Therefore the network is small-world for p < pc while it 
has finite Hausdorff dimension for P > Pc. Moreover the 
maximum degree k^ax increases significantly below the 
transition for P > Pc. 

Furthermore, for P < Pc, the degree distribution P{k) 
is exponential, and the distribution of the curvature P{R) 
has a negative exponential tail,the average curvature is 
{R) = 1/N and its second moment (P^) is finite. 

For p > Pc, instead, P{k) follows a power-law, and P{R) 
has a negative power-law tail. In this case, the average 
curvature is (R) = l/N^ but its second moment (P^) 
diverges. For every value of p, the network has high 
modularity M and a hierarchical structure m with an 
average clustering coefficient C{k) of nodes of degree k 
decaying as C{k) k~^ and a = 1. 

For the Bose-Einstein Network, the most important 
indicator of the phase transition is the maximum degree 
kmax which scales sub-linearly with time for P < Pc and 
linearly for P > Pc, i.e. in this case the most connected 
node is linked to a finite fraction of all the nodes. More¬ 
over, for P < Pc, D increases logarithmically with the 
network size, i.e. the network is small-world, while for 
P > P it decreases significantly. 

Furthermore, for P < Pc, the degree distribution P{k) 
is scale-free, the network has high modularity M, and 
the distribution of the curvature P(P) has a negative 
power-law tail. For P > Pc, instead, P{k) is dominated 
by outlier hubs, the network has low modularity M and 
P(P) has a negative tail dominated by outlier nodes. 

In both cases, the average curvature is (P) = I/A^, and 
its second moment (P^) diverges. The network has a hi¬ 
erarchical structure with an average clustering coef¬ 
ficient C{k) of nodes of degree k decaying as C{k) k~^ 
and a = 1. 


III. QUANTUM NETWORK STATES 
A. The Hilbert space 

Using a similar approach used already in [55H57] . here 
we define quantum network states. In [55H57] an Hilbert 
space is associated to each node and each possible link 
of a network of N nodes. Here we associate to each node 
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FIG. 3: (Color online)Visualization of the Fermi-Dirac Network (with g{uj) given by Eq. (20) and c = 10) for [3 = 0.05,0.5,5 
and N = 1000. For low value of i.e. (3 < fie — 0.14, the network is small-world, for large values of ^d, i. e. {3 > fie — 0.14, 
the network develops a large diameter. The colour indicates the partition into communities found by running the Louvain 
algorithm m- 


p = 0.05 


p = 0.5 


p = 5 



FIG. 4: (Color online) Visualization of the Bose-Einstein Network (with g(uj) given by Eq. (20) and c = 10) for ^ = 0.05, 0.5, 5 
and N = 1000. Eor low value of ^d, i.e. {3 < fie — 0.06, the network is small-world, for large values of /3, i.e. {3 > fie — 0.06, the 
network is condensed and develops a finite diameter. The colour indicates the partition into communities found by running the 
Louvain algorithm m- 


and Hilbert space Hnode and to each link we associate 
two Hilbert spaces l-Liink and l-Lunk- The total Hilbert 
space l-Ljsfoi a network of N nodes is given by 

N N{N-l)/2 N{N-l)/2 

'HN = ^^'Hnode 'Hunk 'Hlink- (21) 

A single realization of a growing geometric network of 
size N in which the nodes are labelled by the time they 
have been added to the network can be mapped to a 
quantum state by mapping the nodes, the links and the 


triangles of the network to quantum states as described 
in the following. 


B. Nodes quantum states 

To every node i we associate a Hilbert space Hnode 
which the Hilbert space of a fermionic oscillator with 
energy uji. Therefore, to every node i of the network we 
associate a node quantum state that can be decomposed 




m= 2 


m= 00 







FIG. 5: (Color online) The maximal distance D from the initial triangle, the maximal degree kmax^ the entropy rate the 
modularity M calculated using the Louvain algorithm and the average clustering coefficient C are plotted as a function of 
the inverse temperature [3 for the Fermi-Dirac Network (m = 2) and for the Bose-Einstein Network (m = oo). The networks 
have nodes with energies following a Poisson distribution g{u}) with average c = 10. The data are reported for networks of 
size N = 10,000 (averaged 30 times), N = 5000 (averaged 60 times) and N = 2500 (averaged 90 times). The predicted phase 
transition for the Fermi-Dirac Network is at /3c — 0.14, for the Bose-Einstein Network it is at /3c — 0.06. 


in the basis 


with Oi = 0,1. The state 




( 22 ) 


\Oi = l,w) 


(23) 
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m = 2 m = 00 







FIG. 6: (Color online) The maximal distance D from the initial triangle, the maximal degree kmax, and the entropy rate 
are plotted as a function of time t for the Fermi-Dirac Network (m — 2) and for the Bose-Einstein Network (m — oo). The 
inverse temperatures are /3 = 0.05 and /3 = 5 respectively, below and above the phase transitions. The networks have nodes 
with energies following a Poisson distribution g(uj) with average c = 10. The data are averaged 20 times. 


is said to contain a particle of energy uj and can be 
mapped to the presence of the node i with energy uj = uji 
in the network. The state 

|oi = 0,cj) (24) 

is an empty state and can be mapped to the absence of 
the node i in the network. In this case the value of uj is 
irrelevant to characterize the state. 


C. Link quantum states 

To every possible link (i,jf) of a network we associate 
an Hilbert space H/in/c-The Hilbert space l-Lunk is chosen 
to be that of a fermonic oscillator. To every pair of nodes 
(i^j) in the network we associate a link quantum state 


that can be decomposed in the basis 

IK)}, (25) 

where Uij = 0,1. The state 

K = 1) (26) 

is said to contain a particle and is mapped to a link (i^j) 
in the network. The state 

\aij = 0) (27) 

is an empty state and is mapped to the absence of a link 
in the network. 
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FIG. 7: Structural and geometrical properties of the Fermi- 
Dirac Network as a function of the inverse temperature [3 
for single network realizations of size N = 10^. The degree 
distribution P{k), the average clustering coefficient C{k) of 
nodes of degree k and the distribution of the curvatures P{R) 
are plotted for (3 = 0.05, 0.5, 5. The networks have nodes with 
energies following a Poisson distribution g(ijj) with average 
c = 10. For low values of i.e. ^ — 0.14, P{k) and P{R) 

are exponential, while for large values of/3, i.e. /3 > /3c — 0.14 
they become power-law. The average clustering coefficient 
C{k) of nodes of degree k always goes like C{k) oc k~^ . 


FIG. 8: Structural and geometrical properties of the Bose- 
Einstein Network as a function of the inverse temperature [3 
for single network realizations of size N = 10^. The degree 
distribution P(k), the average clustering coefficient C(k) of 
nodes of degree k and the distribution of the curvatures P{R) 
are plotted for [3 = 0.05, 0.5,1. The networks have nodes with 
energies following a Poisson distribution g{u) with average 
c = 10. For low values of ^d, i.e. /3 < /3c — 0.06, P{k) and 
P{R) are scale-free, while for large values of ^d, i.e. [3 > 13c — 
0.06, these distributions become dominated by outliers. The 
average clustering coefficient C(k) of nodes of degree k always 
goes like C(k) ex k~^. 


D. Incident triangles quantum state 


To every possible link of a network we associate 
an Hilbert space PLiink- For the Fermi-Dirac network 
state we will assume that this Hilbert space is the one as¬ 
sociated to a fermionic oscillator. For the Bose-Einstein 
network state we will assume that this Hilbert space is 
instead associated with a bosonic oscillator. Therefore to 
each possible link of the network we associate a incident 
triangles quantum state that can be decomposed in the 
basis 


iWij)}, (28) 


with Uij = 0,1 for the Fermi-Dirac Network and Uij = 
0,1, 2,... for the Bose-Einstein Network. The quantum 
number Uij of links for which Uij = 1 , is mapped to 
the number of triangles exceeding one incident to every 
existing link (i, j). Therefore if m = 2 we can only have 
Uij = 0,1 while if m = oo we can have any integer value 
of the occupation number Uij = 0, 1, 2 .... 


E. Quantum network states 

At each time t the quantum network state can be de¬ 
composed into a basis 

|{oi, . (29) 

i i<j i<j 

We consider the following operators that act on the 
nodes, links and incident triangles quantum states. 


F. Creation-Annihilation operators 

The operators are creation-annihilation of 

node quantum states and have anti-commutation rela¬ 
tions 

{6i(w),6](w')} = S{i,j)S{u},u}'), 

{bi{uj),bj{co')} = 0, 

{bl{oj),b]{uj')} = 0, (30) 

where S{x,y) is the Kronecker delta, S{x,y) = 1 for x = y 

and S{x,y) = 0 otherwise. They act on the node states 
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{| 0 i,Wi)} as 


bi{uii) \oi = 0,uii) 
\oi = 1, Wj) 
bl{ui) \oi = 0,Ui) 
bl{uji) \oi = l,u)i) 


= 0, 

= \Oi = 0,U!i) , 
= \Oi = l,U!i) , 
= 0 . 


(31) 


The operators clj,Cij are creation-annihilation opera¬ 
tor of link quantum states and have anti-commutation 
relations 


{ct,ctj = 0. (32) 

where j), (r, s)] = 1 if i = r and j = s or if i = s and 
j = r and (r, s)] = 0 otherwise. They act on the 

link states \aij) as 

Qjf \^ij ~ 0 ) ~ O 5 

^ij \^ij — f) — \^ij — 5 

4j l«ii = 0) = \aij = 1), 

ct \aij = 1) = 0. (33) 


When these operators act on the incident triangles quan¬ 
tum states, they can generate arbitrary occupation num¬ 
bers riij = 0,1, 2,... Their action on the basis {\nij = n)} 
is given by 

dij \nij =n) = ^/n\nij = n - 1), 

=n) = Vn+1 \nij =n + l). 

(37) 


IV. EVOLUTION OF QUANTUM NETWORK 
STATES 

In this section we define a non equilibrium Markovian 
evolution of the quantum network states. The possibility 
of a Markovian evolution of quantum network states is 
not entirely new in the literature, as it has been for exam¬ 
ple proposed in [55] . Therefore we will define a quantum 
network state |'07v(t)) and its evolution with time. The 
quantum network state at every time step can be decom¬ 
posed into the base 

|'0Ar(t)) — ^ ^ ^{oi,uJi,aij,nij} ' 

fof iCLij tTiij } 


Finally we define two classes of creation annihila¬ 
tion operators acting respectively on the incident trian¬ 
gles quantum states of Fermi-Dirac and Bose Einstein 
quantum network states. The creating and annihila¬ 
tion operators d}- and dij acting on incident triangle 
quantum states of Fermi-Dirac network states are anti¬ 
commuting,i.e. 

{dij,4^} = S[{i,j),{r,s)], 

{dij^drpg^ — 0, 

{4,4J = 0. (34) 

When these operators act on the incident triangles 
quantum states, they can only generate occupation 

numbers riij = 0,1. Their action on the basis 

{\nij = 0), \nin = 1)} is given by 

dij \nij =0) =0, 

dij \nij = 1) = \nij = 0), 

dij 4ij = 0) = \nij = 1), 

4 \nij = 1) = 0. (35) 

The creation annihilation operators d\j and dij that are 
acting on the incident triangles quantum states of the 
Bose-Einstein quantum network states have the commu¬ 
tation relations 


d -• d^ 

Liij , 




We start from an initial condition at t = 1 given by 

^ n p4i)44i) 

001,002,003 i=l,2,3 

where Z{1) is fixed by the normalization condition 
('0Ar(l)|'0Ar(l)) = 1- The quantum evolution of the net¬ 
work state is given by a Markov process whose transition 
rate is determined by the unitary operator U 

|'0Ar(t)) = Ut I'lpNit - 1)) , (38) 

where Ut is defined by 


^12^23^13 1^) ’ 


V i,j\i<j 

X&l+2K+2)4+2)iC(t+2)i44^b- (39) 

Here as in Eq. Q Cij = f{uji^ujj) and Z(t) is fixed by the 
normalization condition ('07v(^)|'0Ar(t)) = 1. The quan¬ 
tum operators ht^-, can take two different values, defining 
in this way the Eermi-Dirac Quantum Network state and 
the Bose-Einstein Quantum Network state, i.e. 




d\j Eermi-Dirac Quantum Network state 
d\j Bose-Einstein Quantum Network state 


dij , dj-g 

d^ - d^ 

5 ^rs 


= 0 , 
= 0 . 


In the following section we will consider in detail the 
Eermi-Dirac and the Bose-Einstein Quantum Network 
(3^) states. 
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V. FERMI-DIRAC QUANTUM NETWORK 
EVOLUTION 


A. Path integral 


The evolution of the Fermi-Dirac Quantum Network 
state is given by 

= Ut \tpN{t - 1)) 

V UJt + 2 l,j\l<j 

X&l+2K+2)4+2)iC(t+2)i4cb»J - 1)) (40) 

where eij = and Z = Z{t) is fixed by the nor¬ 

malization condition ('07v(t)|'0Ar(^)) = 1- 

Using the definitions of the creation and annihilation 
operators defined in Sec. |III[ the normalization constant 
Z{t) of the Fermi-Dirac Quantum Network state is fixed 
by the path integral 

= E E (41) 

Wit')} Wt')} 


for t >2. Here is a sequence of links £{t') = 

{it'Wt') and is a sequence of energies of 

the nodes that describes a single history over which the 
path integral is calculated. In the path integral in Eq. 
(41) each path is assigned a weight 


t+2 

= n5(^on 

i=l 


x(l - 


(42) 


where the terms aij{t) and nij{t) that appear in Eq. (42) 
can be expressed in terms of the history {i{t')}t'<t as 


aij (t) 




t-i 


E + 2,it')] +S[{i,j),{t' + 2,it')]) 

t'=2 


+6[{i,j), (1,2)] + (1,3)] + (2,3)], 


t-1 


E-^[^(^'),(i,i)]- 

t'=2 


(43) 


Therefore Z{t) can be interpreted as a partition function 
of a statistical mechanics problem in which each path up 
to time t has probability 




m 


(44) 


Each of the paths {uj(t')^£{t')}t'<t can be mapped to 
a geometrical network evolution with m = 2. In this 
mapping uj{t) indicates the energy of the node added to 
the network at time t, £{t) = {itWt) indicates the link to 
which we attach a new triangle at time t, aij{t) indicates 


the adjacency matrix of the network, nij{t) indicates the 
additional number of triangles incident to an existing link 
(i, jf). The probability of each geometrical network evo¬ 
lution described in Sec. IIA is the 

same as the weight that the corresponding history has in 
Eq. 

Starting from Eq. (44) we can calculate the condi¬ 
tional probability that at time t we add a link £{t) = 
given the present state of the network evolution, 

n[tj) = A straightfor- 

ward calculation shows that 


n 


[F] 

dd) 


(t) 


Zf 


(45) 


where 


ZF = Ze - «ii(*))- (46) 

i<j 

Given that 1 — riij = 1 has the graphical network inter¬ 
pretation ^ij = 1, i.e. indicates that the link (i, j) is not 
saturated, while I — riij = 0 indicates that the link is 
saturated, i.e. ^ij = 0, the expression in Eq. ( [4^ is the 

same as defined in Eq. (j^. It follows that studying 

the geometrical network evolution for m = 2 determines 
the properties of the Eermi-Dirac quantum network state. 
Eor this reason we call the growing geometrical network 
with m = 2 the Eermi-Dirac Network. 


B. Fermi-Dirac Statistics 

The average of the quantum number riij over all the 
links of the Fermi-Dirac Network follows the Fermi-Dirac 
distribution. Since riij = O 5 1? equivalently, we can say 
that the probability that a link with energy e is saturated 
follows the Fermi-Dirac statistics. To derive this result, 
let us consider the master equation [4] for the number 
Np{n\ujw') of links (i, j) (with i > j^ uji = uj and ujj = 
uj')^ that have riij = n = 0,1 at time t. Since at each time 

we choose a link (i, j) with probability (t), only if 
it is unsaturated (i.e. riij = 0 ), we add one triangle to 
the link (i.e. riij = 0 ^ riij = 1 ) and we add other two 
unsaturated links, the master equation reads 

p-/3e 

NP~^{n = l\ujw') = — Np{n = 0\ujw') 

Zf 

-\-Nf{tI = l\idw') 

g-/3e 

NP~^{n = 0\uj,uj') = —— NF{n = 0\uj,uj') 

Zf 

+2pf{oj, uj') + NF{n = 0|u;, cj'), (47) 

where e = and pf{ojw') Is the probability that 

a new link (i,j) of the network with i > j links two 
nodes with energy uji = uu and ujj = uj'. In order to solve 
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this master equation we assume that the normalization 
constant Zp oct and we put 


with probability e it follows that the recursive 

equation for gri^) reads 




g-/ 3 MF = lijjj 

t—yoc t 


(48) 


This is a self-consistent assumption that must be verified 
by the solution of Eqs. ( [47| ). Moreover we also assume 
that at large times Np{n\LJ^uj') ^ 2tPF{n\uj^uj'). In fact 
the number of links in the network is 2t + 1 2t for t ^ 
1. Here Pi?(n|u;, cj') indicates the asymptotic probability 
that a random link (i,j) with i > j and uJi = uj, ujj = uo' 
has riij = n. With these assumptions, we can solve Eqs. 
(47) finding 


Pp(n = Olcj, cj') = 

= pf{uj,uj') 
PF{n = = pf{oj^uj') 


= pf{uj,uj')[ 1 - npie)] 

1 

g/3(e-AXF) 1 
= pF{^,^')nF{^) 


(49) 


where e = and ni?(e) is the Eermi-Dirac occupa¬ 

tion number 

ufU) = ^7—^(50) 

Considering all the links with energy e, we have 

(n|e) = pF{e)nF{e), (51) 

where 

pFie) = '^S[e,f{u),uj')]pFi<^,u)') (52) 

UJ,UJ' 

and where 


PFico) = [Pi.(n = 0|a;,w') 

UJ' 

+Pp(n = 0|cj',u;)] 

= +PF(^',^)]riF[/(^,^')](55) 


where in the last equation we have used the expression 
for Pf{ti = Olci;,^^') given by Eq. (49). Eq. (55) can be 
formulated as the eigenvalue problem 


5f(w) = y] Af(w,w')5(w'), 


(56) 


where 


Af{lo,u') = g{u)nF[f{to,u')] 

X jl -'^g{uj")nF[f{oJ,uj")] 


-1 


(57) 


Since we require that gF{^) is a probability, i.e. it is non¬ 
negative and normalized, the solution of the eigenvalue 
problem is given by the Perron-Erobenious eigenvector 
gF{^) of the matrix Af{oo,uj') satisfying 


Z,9f{uj) = 1. 


(58) 


Einally the chemical potential /ii? is fixed by the self- 
consistent condition in Eq. (48) that can be rewritten 
as 


'Y^PF{e)nF(e) = - 


(59) 


{n\^) = Z,S[e,f{LO,uj')]pF{uJ,uj') nPF(n|w,w'). (53) 

uj,uj' n=0,l 

Therefore, in the Eermi-Dirac Network the average of the 
incident triangles quantum number over links of energy 
e follows the Eermi-Dirac distribution. 

To complete the solution it is necessary to find the 
correct expression for pf{oo^uj'). Since by definition uj is 
the energy of the new node attached to the network at 
time t, and since this energy is drawn randomly from a 
distribution g{uj)^ we have that the probability pf{oj^uj') 
can be factorized, 

/)f(w,w') = fl(w)flF(w'), (54) 

where gF{^') is the probability that a new triangle is at¬ 
tached to a link having at its end a node of energy uj' and 
is thus normalized. Therefore we can write a recursive 
equation for gF{oj'). Since we attach new triangles to 
random unsaturated links with energy eij = e = f{uj,uj') 


which is the same equation as the one fixing the chemical 
potential in a Eermi gas m with density of states Pf{^)i 
inverse temperature p and specific volume v = 2. If the 
self-consistent equation given by Eq. (48) has a solution, 
and Zf oc the master equation asymptotically in time 
has a stationary solution given by Eqs. (49). This implies 
that the average of the quantum numbers riij over links of 
energy e, i.e. (n|e), follows the Eermi-Dirac distribution. 


C. Structural properties of the Fermi-Dirac model 

Let us here characterize some of the important struc¬ 
tural properties of the Eermi-Dirac network model. Eirst 
of all, let us consider the degree distribution P{k). In 
order to find P{k) we first write the master equation for 
the number Njp{k\uj) of nodes that at time t have degree 
k given that they have energy uji = uj. Eor simplicity 
in this paragraph we consider the linear relation Eq. § 
between link and node energies. 
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The master equation [4] for Np{k\uj) reads 


N^+\k\uj) = - - - N^k - l\uj)[l - 5{k,2)] 

p-Piuj-fip) 

- - - Np{k\uj) + g{io)5{k, 2) 

+N],{k\w) (60) 


where we have assumed that asymptotically in time we 
can define the chemical potential flp given by 


g/3/iF = ijjjj 

t—yoo 


'J2ije l^'^^{l-nij)aijS{ki,k) 

J2iHki,k) 


By assuming in the large network limit t ^ 1 that 
NF{k\uj) ^ tP{k\uj)^ solving Eq. (60) we get, 


P{k\uj) = g{L 0 ) 




+ l] 


k-l 


(61) 


for k > 2. Therefore, summing over all the values of the 
energy of the nodes uj we get the full degree distribution 
p{k) 


P{k) = Y,g{u:) 

UJ 


^^{uj-flF) 

\^^l3{uj-jlF) + ij ^ ^ 


(62) 


for k >2. 

The curvature Ri of a node i is given by 


Ri = 


4 - 
6 


(63) 


So the distribution of the curvature P{R) is given by 


P(i?) = ^5H 


^^/3{uJ-flF) l] 


3(l-2it:) 


(64) 


where R < Therefore the distribution of the curva¬ 
ture is decaying exponentially for negative values of the 
curvature.Moreover the average curvature is (R) = 1/N 
and the fluctuations around this average are bounded,i.e. 
(R^) < oc. 


D. Comparison with numerical simulations 


Here we numerically simulate a Fermi-Dirac Network 
Evolution in which the energies of the nodes are non¬ 
negative integers with g{uj) given by a Poisson distribu¬ 
tion with average c as in Eq. ( [^ and link and node 
energies related linearly Eq. <§• 

We compare the results of the theory with the out¬ 
comes of the simulations as long as the chemical potential 
fiF defined in Eq. (59) is well defined. In particular we 


show the results of simulations confirming that the aver¬ 
age of the quantum number riij over the links of energy 
e, {n\e) is given by Eq. (51) and follows the Fermi-Dirac 


statistics. In Figure we compare the results of the sim¬ 
ulations with the theoretical expectations by plotting the 
right and left hand side of equation 


nF{e) 


{n\e) 

PF{e) 


(65) 


equivalent to Eq. (51) and the degree distribution P{k) 
of this network with the theoretical expectation given by 
Eq. ( [^ . We find very good agreement in both cases as 
displayed in Figure pT) 


E. Phase transition in the Fermi-Dirac network 


The self-consistent approach for solving the Fermi- 
Dirac network is based on the assumption that the chemi¬ 
cal potential /ii? of the network defined in Eq. (48) exists. 
But in the network it is possible to find a phase transition 
at high enough inverse temperat ure, i .e. for P > Pc where 
this assumption fails (see Sec. IIC). In order to deter¬ 
mine where this phase transition occurs we have solved 
the self-consistent equation for the chemical potential juf 
given by Eq. (59). This equation can always be solved 
to find the chemical potential /ii?, but the value of the 
chemical potential jif as a function of the inverse tem¬ 
perature can have a maximum for P = Pc- Here we have 
identified this maximum with the onset of the phase tran¬ 
sition. In fact, from the dynamic rules of the model it is 
clear that the network dynamics for increasing value of P 
tends to attach new triangles on unsaturated nodes with 
lower energy. Therefore if the probability P{n = Olcc;,^^') 
is given by Eq. (49) the chemical potential fiF that can 
only increase with increasing inverse temperature p. 

In Figure [TT] we show the chemical potential /ii? as a 
function of the inverse temperature p for a Fermi-Dirac 
network with a Poisson distribution g{ijj) (given by (20)) 
with average c = 10 and linear relation between the en¬ 
ergy of the nodes and the energy of the links (Eq. ([^). In 
order to perform the numerical calculation of the chem¬ 
ical potential jiF the distribution g{uj) is truncated at a 
cutoff value cja = 100. The chemical potential ^f has 
a maximum at Pc — 0.14 which is a good prediction for 
the phase transition as can be seen from the simulation 
results shown in Figure 


VI. BOSE-EINSTEIN NETWORK EVOLUTION 
A. Path integral 


The evolution of the Bose-Einstein quantum state is 
described by the unitary operator Ut defined in the fol- 
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FIG. 9: (Color online) The Fermi-Dirac occupation number nrie) is extracted from the simulation results by plotting {n\e)/prie) 
(star points) and compared with the theoretical prediction (solid red line). Data are shown for a Fermi-Dirac network with 
g{(jo) given by a Poisson distribution with c = 2. The network size is N = 2000 and the simulation results are averaged over 30 
runs. 
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FIG. 10: (Color online) The degree distribution P{k) of the Fermi-Dirac Network is plotted with star points for a Fermi-Dirac 
network with g{(jo) given by a Poisson distribution with c = 2. The network size is N = 5000 and the s imu lation results are 
single realizations of the networks. The simulations are compared with the theoretical predictions of Eq. (62) shown as a solid 
red line. 
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FIG. 11: (Color online) The chemical potential fip versus 
the inverse temperature [3 for a Fermi-Dirac Network with 
g{uj) = with (jj G [0,100] , c = 10 and with JV indicat¬ 

ing the normalization sum. The critical value of the inverse 
temperature is /3c — 0.14 which is in good agreement with the 
simulations (see Figure 0). 



lowing, 


|'0Ar(^)) = Ut |'0Ar(^ - 1)) 

V <^t+2 i.,j\i<j 

>^K+2i^t+2)c\^^2)At+2)Aj4j(^ij - 1 )) ,( 66 ) 

where €ij = and Z = Z(t) is fixed by the nor¬ 

malization condition ((ijj nn{ t)) = 1- 

In this case the normalization constant Z{t) is given 
by the path integral 

= E E (67) 


for t > 2 and = f{uji,ujj). Here 
is a sequence of links £{t') = and a sequence of 

energies of the new nodes uj{t') that describes a single 
history over which the path integral is calculated. Each 
path {uj(t')^£{t')}t'<t ia Eq. (67) is assigned a weight 


x(l + n^((,)(i'))e 


t+2 t 

) = n 

2 = 1 t'=2 


( 68 ) 


where aij(t) and nij(t) can be expressed in terms of the 
history {£{t')}t'<t in the same way as in Eq. (43). Note 


the characteristic sign difference in Eq. ( |68[ ) compared to 
the Eermi-Dirac case. Therefore Z{t) can be interpreted 
as a partition function of a statistical mechanics problem 


in which each path up to time t has probability 




m 


(69) 


Each of the paths {uj{t')^£{t')}t'<t can be mapped to 
a geometrical network evolution with m = oo. In this 
mapping uj{t) indicates the energy of the node added to 
the network at time t, i{t) = indicates the link 

to which we attach a new triangle at time t, and aij{t) 
indicates the adjacency matrix of the network. Moreover 
nij{t) indicates the number of triangles that exceed one, 
attached to the link (i, j). The probability of each geo¬ 
metrical network evolution described in 

Sec. |II A| is the same as the probability given by in Eq. 


Starting from Eq. (68) we can calculate the condi¬ 


tional probability that at time t we add a link £{t) = 
{i^j) given the present state of the network evolution, 

A straightfor- 

ward calculation shows that 


n!3)(^) = 


e 






(70) 


where 


ZB = J2e-^^- 

i<j 


aij{t){l + riijit)). 


(71) 


The expression in Eq. (45) is the same as n[G.., defined 
in Eq. ^ for m = oc. In this case = 1 for nodes 

{i,j) for which there is a link, i.e. ciij{t) = 1. It fol¬ 
lows that studying the geometrical network evolution for 
m = oc determines the properties of the Bose-Einstein 
quantum network state. Eor this reason we call the grow¬ 
ing geometrical network with m = oo the Bose-Einstein 
Network. 


B. Bose-Einstein statistics 


In the Bose-Einstein Network the average of the quan¬ 
tum number riij over links with energy e follows the Bose- 
Einstein distribution. To show this result let us consider 
the master equation [4] for the number of 

links (i^j) (with i > j and Ui = u and ujj = uj') that 
have riij = n = 0,1, 2 ,... at time t. Since at each time 

we choose a link {i,j) with probability the master 

equation reads 


Vy^(n|a;,a;') 




D-/3e 


-N^{n — l|cc;,cc;')[l — (^(n,0)] 


Zb 


-\-2pB{co,ui')S{n,0) Ng(n\ui,ui'), 

(72) 
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where e = and pb{oj^uj') is the probability that 

a new link will connect a new node with energy uj and 
an old node with energy uj'. In order to solve this mas¬ 
ter equation we assume that the normalization constant 
Zb oc t and we put 


e 


t—yoc t 


(73) 


Moreover we also assume that at large times 
Nh {n\uj,uj') 2tPB{n\uj,uj') as the number of links in 
the network is 2t + 1 2t for t ^ 1. The quantity 
Pb{ti\uj,uj') indicates the asymptotic probability that a 
random link (i^j) with i > j and uJi = uj and ujj = uj', 
has riij = n. Making these assumptions, it is possible to 
solve Eq. ( [7^ getting 


PB(n, 


.. r(n + l) 
r [e/^h-Ms) + n + 2] 


for n > 0. Eq. ( [74| ) is automatically normalized once the 
distribution pb(u^uj') is normalized. Therefore the aver¬ 
age of the quantum numbers riij over links with energy 
eij = e is given by 

{n\u,Lo') = 

= />B(w,w')«B(e), (75) 


where riBi^) indicates the Bose-Einstein occupation num¬ 
ber 


nsie) 


1 

— I ' 


(76) 


Summing over all the links with energy e, we get 


(n|e) = pB{e)nB{e), (77) 


where 


PB{e) = '^5[e,f{Lo,LoJ]pB{uj,uj'). (78) 

UJ,UJ' 


Therefore in the Bose-Einstein network the average of 
the quantum number riij over links of energy e follows 
the Bose-Einstein distribution. 

To complete the solution it is necessary to find the cor¬ 
rect expression for pb{uJ:Uj')- Since by definition uj is the 
energy of the new node attached to the network at time 
t, and since this energy is drawn randomly from a distri¬ 
bution g{uj), the probability pb{jo^Jo') can be factorized, 

Pb{uj,uj') = g{uj)gB{uj'), (79) 

where gB{^') is the probability that a new triangle is 
attached to a link having at its end a node of energy uj' 
and is therefore normalized. 


We can write a recursive equation for gB{uj')- In fact 
we have 

gB{oj) = {[pb(w,w') + (n|w,w')] 

Uj' 

P[Pb{(^' + {n\uj'^uj)]} 

= y][pB(w,w') + PB(w',w)]nB[/(w,w')], (80) 


where in the last equation we have substituted Eq. (75) 
into {n\uj,uj'). This equation can be formulated as the 
eigenvalue problem 


5b(w) = y] Ab(w,w')5(w'), 


(81) 


where 

AB{oJ,ujj = g{oj)nB[f{uJ,ujj] 


-1 


X 11 - E [f{0J + w")] I . (82) 

Since we require that gB{^) is a probability, i.e. it is non¬ 
negative and normalized, the solution of the eigenvalue 
problem is given by the Perron-Erobenious eigenvector 
gB{^) of the matrix Ab{uj^uj') satisfying 


E5b(w) = 1. 


(83) 


By imposing the self-consistent condition in Eq. (73) we 
find the equation determining the chemical potential p b , 

Y.PB{e)nB{e) = \, (84) 

e 

which is the same equation as the one fixing the chemical 
potential in a Bose gas m with density of states Pb{^)i 
inverse temperature /3 and specific volume v = 1/2. 

If the self-consistent Eq. (86) has a solution, this im¬ 
plies that the average of the quantum number riij over 
links of energy e, (n|e), follows the Bose-Einstein distri¬ 
bution. 


C. Structural properties of the Bose-Einstein 
Network 

In this section we derive the degree distribution and 
the distribution of the curvature in the case in which 
the energies of the links are linearly dependent on the 
energies of the nodes as in Eq. ^ . In order to derive the 
degree distribution in this case, let us write the master 
equation [4] for the number N'^{k\uj) of nodes that at 
time t have degree k and energy cj, i.e. 

^ - >-Npk-l\co)[l-6{k,2)] 

p-P{u)-fLB)p 

-^- Npk\uj) + gUo)5{k, 2) 

+N%{k\L0), 


(85) 
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where we have assumed self-consistently that asymptot¬ 
ically in time /Ib is defined as 

( 86 ) 

Assuming that asymptotically in time N'^{k\uj) 
tP{k\uj) we find the expression for P{k\uj) and substi¬ 
tuting this expression in Eq. (85) we have, 


^ T + l] ^ ' 

for k > 2. Therefore the degree distribution of the entire 
network is scale-free and given by 


P(fc) = ^P{k\co) 


( 88 ) 


for k > 2. The fraction of the total number of links 
attached to nodes with energy uj is given, asymptotically 
in time, by 

(^1^^) = ^ E (89) 

k 


where 7 x 3 ( 00 ) is defined as 

nsioj) = 


1 


— 1 ' 


(90) 


In Eq. (89) the first term indicates the fraction of links 


initially attached to the new nodes of energy cj, and is 
therefore given by g(uj) because every new link has one 
end attached to a new node and the new node has energy 
00 with probability g{oo). The second term represents the 
fraction of links attached to the nodes of energy cj, after 
the time of their arrival into the network. This term is 
proportional to the Bose-Einstein occupation number. 

The curvature Ri of a node i is given by 


Ri = 


4 - 
6 


(91) 


So the distribution of the curvature P(R) is given by 

^ y. fl(c8;)r[2 + e/^(--AB)]e/3(u7-AB)r(4 - 6P) 

^ ^ T + 5 _ ^ ^ 

where R < ^. Therefore the distribution of the curvature 
is scale-free and decaying as a power-law for negative 
values of the curvature.Moreover the average curvature 
is (R) = 1/N and the fluctuations around this average 
N ^ oc are diverging, i.e. {R‘^) ^ 00 as ^ 00 . 


D. Comparison with numerical simulations 

We numerically simulate a Bose-Einstein Network Evo¬ 
lution in which the energies of the nodes are non-negative 


integers with g(oo) given by a Poisson distribution with 
average c (Eq. (2^) and the link and node energies are 
related by Eq. ( ^ and compare with the theoretical re¬ 
sults. We verify that Eq. (77) is satisfied by the results 


of the numerical simulations. In Eigure[T^ we plot both 
sides of the equation 


nsiuj) = 




(93) 


equivalent to Eq. (7^ and the degree distribution P{k) 
of this network with the theoretical expectation given by 
Eq. ( 88 ), finding very good agreement between theoret¬ 


ical and numerical results in both cases (see Eigure 13). 


E. Transition: Bose-Einstein condensation 

The self-consistent approach for solving the Bose- 
Einstein network is based on the assumption that the 
chemical potential /tb of the network defined in Eq. (73) 


exists. But in the network it is possible to find a phase 
transition at high enough inverse temperature, i.e. for 
P > Pci where this assumption fails. Since the negative 
chemical potential, /iB < 0 , can only increase with the 
temperature, the critical value of the inverse temperature 
Pc is determined by the self-consistent equation for the 
chemical potential Eq. (84) where we impose /iB = 0 and 
p = Pc, i.e. 




1 


Pee _ 1 


(94) 


Here p%(e) is given by Eq.(79) and is calculated for /iB = 
0 . In a Bose Gas of density of states Pb{^) the existence 
of a finite critical temperature Pc indicates the onset of 
the Bose-Einstein condensation. Eor the network this 
indicates that there is a single link incident to a finite 
fraction of all the triangles, and therefore also the degree 
of the incident nodes is a finite fraction of the total degree 
of the network (see Eig ures ( ^, ([^, ^ and the discussion 
of the transition in Sec. IIC). This phenomenon is similar 


to the one occurring in other models displaying the so- 
called Bose-Einstein condensation in complex networks 
min]. In Eigure we show the chemical potential ps 
as a function of the inverse temperature P for a Bose- 
Einstein network with a Poisson distribution g(uo) with 
average c = 10 (given by Eq. ([^) and a linear relation 
between the energy of the links and the nodes (given 
by Eq. (§), where in order to perform the numerical 
calculation of /iB the distribution g(uo) is truncated at a 
cutoff value uo\ = 100. In this network a Bose-Einstein 
phase transition occurs at /dc ^ 0.06 which is in very good 
agreement with the simulation results shown in Eigure 
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FIG. 12: (Color online) The Bose-Einstein occupation number ns(e) is extracted from the simulation results by plotting 
(n|e)/ps(e) (star points) and compared with the theoretical prediction (solid red line). Data are shown for a Fermi-Dirac 
network with g{oj) given by a Poisson distribution with c = 2. The network size is iV = 1000 and the simulation results are 
averaged over 50 runs. 




FIG. 13: (Color online) The degree distribution P{k) of the Bose-Einstein Network is plotted with star points for a Bose- 
Einstein Network with g{oj) given by a Poisson distribution with c = 2. The network size \s N — 5000 and the simulation 
results correspond to a single network realization. The simulations are compared with the theoretical predictions of Eq. (88) 
shown as a solid red line. 


VII. THERMODYNAMICS OF QUANTUM 
GEOMETRIC NETWORKS 


and the entropy S{t) of the quantum geometric network 
evolution as 


A. Relation between the total energy E and the 
entropy S 

Given the quantum geometric network evolution it 
is natural to characterize its thermodynamic properties 
above and below the phase transition. Let us define the 
total energy of a quantum geometric network as 

E{t) ( 95 ) 

i<j 


s{t) = - E 

t' <t 

(96) 


In this expression, is the prob¬ 

ability that the temporal evolution of the network until 
time t is described by the subsequent addition of trian¬ 
gles to the links given that the energies of 

the nodes until time t — 1 are Together with 

the definition of given by Eq. (11) it can be easily 
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FIG. 14: (Color online) The chemical potential /ib versus 
the inverse temperature 13 for a Bose-Einstein Network with 
g{uj) = with CjJ G [0,100], c = 10 and Af indicating 

a normalization sum. The critical value of the inverse tem¬ 
perature is ydc — 0.06 which is in good agreement with the 
simulations (see Figure]^. 


derived that 

if[1](f) = AS{t) = S{t) - S{t - 1). (97) 

Moreover we have already found (Eq. ( |12| )) that 

if'^'=/3(ei,)+In Z, (98) 


where the average (e^j) is given by Eq. (13) and can be 
related to the expected increment in time of the energy 
E{t) given by Eq. ( [^ , 

(e^,-) = (AE) = {E{t) - E(t - 1)). (99) 

The relation between AS' and {/AE) can be found using 
Eqs. (97) — (98) and (99). This relation depends on the 
inverse temperature For /3 < /3c we have that for large 
times, t ^ 1, Z and therefore 


AS P [{/AE) — /i] + Int. 


( 100 ) 


In the limit t ^ oc, the chemical potential converges 
to jj^F for the Fermi-Dirac network and to for the 
Bose-Einstein network. Instead, for P > Pc we have that 
Z = 0(1). By putting Z = = 0(1), we have 

AS = P[{AE)-u] (101) 

where z/ is a stochastic variable depending on the history 
of the network. 


B. Probability of an unlabeled quantum network 
final state 

Given a network quantum state mapped to a geometric 
network with A = t + 2 nodes i = 1, 2 ,..., A, 

= 11 \L 0 i) 11 \aij) 11 \rnj) , ( 102 ) 

i i<3 i<j\aij = l 


consider the unlabelled quantum network state con¬ 
structed from it by considering all the permutations of 
the node labels 


\ps{'/)) — rr rr 

{7r(i)} i i<j 




(O' 


r0)>, (103) 


X n 

*<i|a7r(07r(i)=l 

where 7r(i) indicates a permutation of all the indices 
{i = 1, 2 ..., A} of the nodes. We are interested in evalu¬ 
ating the probabilit y Pf that the quantum network state 
|'07v(^)) given Eq. (38) is found in this unlabelled final 
network state, i.e. we want to calculate 

Pf = \{<i>smNm^- ( 104 ) 


This is equivalent to calculate the probability that the ge¬ 
ometrical network evolution generate a time t networks 
that are equivalent under relabeling to the nodes. Any 
given history from time t = 1 to time t, corresponds to 
a series of processes consisting in gluing new triangles to 
unsaturated links. Therefore, given a final network state 
one can prune the network in the same order at which the 
triangles have been added. In this pruning process, we 
start by removing the last triangle it has been glued to 
the network, then removing the second last triangle un¬ 
til we reach the initial triangle of the network evolution. 
Nevertheless, different network evolutions can generate 
networks that only differ by a relabeling of the nodes 
corresponding to the same set of triangles added in a dif¬ 
ferent order starting from the initial triangle. Finding 
how many such histories exist is a problem can be cast 
into the problem of finding the number N of ways we can 
prune the network corresponding to the final unlabelled 
network state. The pruning of a given geometric net¬ 
work consists in an iterative process in which one takes 
randomly any triangle incident only to a single other tri¬ 
angle and different from the initial triangle, and removes 
it, until the full network is reduced to the single initial 
triangle. If we call Af the number of ways this pruning 
can be done, we find that 


Pf = 




(105) 


where E is given by Eq. ( [^ . Finding A is a combina¬ 
torial problem to be solved for any given final network 
realization. The formula for Af can be found iteratively 
using similar techniques as in uniiHi]- We can interpret 


>S = 21og(A/') 


(106) 


as an entropy associated to the unlabelled network state. 
With this notation we have for the Fermi-Dirac Network, 

Pf = (107) 

and for the Bose-Einstein Network 

Pf = -:z-, 


(108) 
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VIII. GENERALIZED GEOMETRIC NETWORK 
MODEL WITH ENERGY OF THE LINKS 

A. The network evolution 


Here we consider an extension of the geometric network 
model that might allow the generation of networks in 
which not all the nodes are at the boundary. In particular 
our goal is to describe geometric network models in which 
the resulting networks might contain saturated nodes, i.e. 
nodes incident only to saturated links. 

Therefore here we define a generalized geometric net¬ 
work model where the non equilibrium network evolution 
includes and additional processes with respect to the one 
introduced for the geometric networks studied in the pre¬ 
vious sections. We start from a network formed by a sin¬ 
gle triangle, a simplex of dimension = 2. Each link 
can belong at most to m triangles. If a link (i, j) belongs 
to m triangles it is saturated and ^ij = 0. If it belongs 
to less than m triangles it is unsaturated and ^ij = 1. 
To each node i an energy of the node uJi is assigned from 
a distribution g{uj). The energy of the node is quenched 
and does not change during the evolution of the network. 
Moreover to each link (i, j) we associate an energy of the 
link eij given by a symmetric function of the energy of 
the nodes i and j as in Eq. 0. 

At each time we perform two processes: process (a) 
and process (b). The process (a) is the same process 
considered in the original model described in Sec. ra 
Here we consider also an additional process (process (b)) 
occurring at each time with probability p. Process (a) 
and process (b) are described in the following. 

• Process (a)- We add a triangle to an unsaturated 
link (i, j) of the network linking node i to node j. 

We choose this link with probability given 

by Eq. Having chosen the link (i,j) we add 

a node r, two links (i,r) and (j, r) and the new 
triangle linking node i, node j and node r. 


• Process (b)- With probability p we add a single 
link between two nodes not already linked and at 
hopping distance 2, and we add all the triangles 
that this link closes, without adding more than m 
triangles to each link. In order to do this, we define 
a variable cr(t) = 1 if process (b) takes place at 
time t (event which occurs with probability p) and 
cr(t) = 0 if process (b) does not take place at time 
t. If cr(t) = 1 we choose two unsaturated links 
(i',/) and (/,r') specified by q = with 

probability given by 
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FIG. 15: (Color online) The generalized growing geometric 
network model is the underlying network of a growing sim¬ 
plicial complex evolving by process (a) and process (b). In 
process (a) a new triangle is connected to the network and 
glued to an existing unsaturated link of the network. In pro¬ 
cess (b) two nodes at distance 2 connected by unsaturated 
links are connected by a new link, and all the triangles that 
this link closes are added, provided that no more than m tri¬ 
angles are incident to each link. In the figure the case with 
m = 2 is plotted. 


where C is the normalization constant and where 
^ij = 1 indicates an unsaturated link, while ^j j = 0 
indicates a saturated link. Moreover in Eq. (109) 
the quantity 1 + riij indicates the total number of 
triangles incident to the link (i, j). Then we add a 
link (i', r') and all the triangles that the link (i', r') 
closes. 


In Eigure 15 we show how the generalized growing net¬ 
work model can be extracted by the generalized growing 
simplicial complex evolving through process (a) and pro¬ 
cess (b). 

In Eignrep!6|we show schematically the dynamical rules 
for building generalized geometrical growing networks 
with m = 2 and m = oc that we will call respectively 
Generalized Eermi-Dirac Network and Generalized Bose- 
Einstein Network. 


B. Entropy rate of the Generalized Geometric 
Network 

Calling = {ci;(t),-^(t), cr(t), q(t)} the generalized 

geometric network evolution is described by the sequence 
At time t the sequence includes the new 
terms Q{t). The entropy rate of the generalized geometric 
growing network is 

Hoit) = 

n(t) q(t) 


X 


( 110 ) 
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FIG. 16: The Generalized Fermi-Dirac Network evolution 
(panel A) and the Generalized Bose-Einstein Network Evolu¬ 
tion (panel B). At each time a new triangle is added to a link 
{i,j) chosen according to the probability given by Eq. 

(process (a)), and with probability p also process (b) takes 
place and a new link is added between two nodes at dis tance 
2 with probability chosen with given by Eq. (109). The 
maximal number of triangles incident to a link are m — 2 for 
the Generalized Eermi-Dirac Network Evolution and m — oo 
for the Generalized Bose-Einstein Network evolution. In the 
figure, at time t = 3 the process {h) takes place and the new 
link added according to this process is plotted with a thick 
green line. 


where 7^(f](t)|{f](t')}t'<t)) indicates the probability of 
having 0(t) given {0(t')}t'<t. The probability of uj{t) = 
uj is and it is independent of all the other precedent 
events. The probability of i{t) is given by as in Eq. 

The probability of having process (b), i.e. a{t) = 1, 
is p, and the probability of not having process (b), i.e. 
a{t) = 0, is 1 — p. Finally if the process (b) takes place, 
the probability of q(t) is as in Eq. (109). The 

entropy rate of the network evolution can therefore be 
written as the sum of three different entropy rates 


Hcit) = 


(Ill) 


Here, is the contribution to the entropy rate due to 
the time-independent random distribution of the energy 


of the nodes, 

( 112 ) 

UJ 

The quantity indicates the contribution due to 

process (a). 


i/W(i) = -^nW^innW^, 

m 


(113) 


while iJPl(t) indicates the contribution due to process 

(b), 


ij[2l(i) = -(1 - p) ln(l - p) - 


q(^) 


P%it) 


.(114) 


A change in the scaling of (t) + (t) with time in¬ 
dicates a phase transition in the network. Such a phase 
transition can occur at high values of the inverse tem¬ 
perature /3, where the network dynamics can become ex¬ 
tremal, similar to what we have seen occurring in the 
precedent sections in the case p = 0. 


C. Phase transition in the generalized geometric 
networks 

Except for the case m = 2 which is planar, the gener¬ 
alized geometric network model with p > 0 is not planar, 
and we have 


lim - > 0. (115) 

t^oo t 

In fact, if process (6) occurs, we add zero new nodes, one 
new link, and we have the possibility to add more than 
one triangle if m > 2. Therefore we have that at any 
given time 


(Ax) = (x(i)-x(i-1)) > 0. (116) 


For these networks we extend the definition of the cur¬ 
vature defined for planar graphs [ISIIIT] and we take the 
curvature Ri associated to each node i of the network 
given by 

= + ( 117 ) 


which is the same definition as in Eq. (19) apart from 
that there is no simple relation between ki and Ti. Here 
we focus in particular on the cases m = 2 and m = oc 
and we show numerical evidence that a phase transition 
might occur in these networks. Specifically, we consider 
the case in which uj can only take integer values and the 
distribution g{uj) is Poisson with average c, i.e. Eq. ( [^ , 
and node and link energies are related by Eq. (§. 

As a function of (3 we observe a phase transition both 
in the case of the Generalized Fermi-Dirac Network and 
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FIG. 17: (Color online) Visualization of the Generalized Fermi-Dirac Network with a Poisson energy distribution of the nodes 
g{uj) given by Eq. (20), c = 10, /3 = 0.05,0.5,5, N = 1000 and p = 0.9. For low value of /3, i.e. [3 < the network is 
small-world, for large values of /3, i.e. /3 > /3c, the network develops a large diameter. The color indicates the partition into 
communities found by running the Louvain algorithm iza. 



FIG. 18: (Color online) Visualization of the Generalized Bose-Einstein Network with a Poisson energy distribution of the nodes 
g(uj) given by Eq. (20), c = 10, yd = 0.05, 0.5, 5,V = 1000 and p = 0.9. Eor low value of yd, i.e. yd < the network is 
small-world, for large values of yd, i.e. yd > I3c, the network is condensed, and develops a small diameter. The colour indicates 
the partition into communities found by running the Louvain algorithm Eg. 


in the case of the Generalized Bose-Einstein Network. 
For f3 > f3c the structure of the network and its ge¬ 
ometry change drastically as it can been seen already 
from the visualizations of the networks (Figure pT] for 
the Generalized Fermi-Dirac Network and Figure |T^ for 
the Generalized Bose-Einstein Network). The transi¬ 
tion is characterized by a different scaling of the en¬ 
tropy rate below and above the tran¬ 
sition. For yd < ydc, increases with 

time as ln(t). For jS > pc instead, 

= 0(1) fluctuates widely during the 


network evolution. Here we discuss in detail the conse¬ 
quences of this phase transition in the Generalized Fermi- 
Dirac Network and in the Generalized Bose-Einstein Net¬ 
work. In Figurewe show major geometrical and struc¬ 
tural properties of the network as a function of the inverse 
temperature yd across the phase transitions. In particular 
we display the maximal shortest (hopping) distance from 
the initial triangle D, the maximum degree kmax of the 
network, the entropy rate the modularity M 

calculated using the Louvain algorithm Eg and the aver¬ 
age clustering coefficient C across the phase transitions. 
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FIG. 19: (Color online)The maximal distance D from the initial triangle, the maximal degree kmax, the entropy rate 
the modularity M calculated using the Louvain algorithm and the average clustering coefficient C, are plotted as a function 
of the inverse temperature /3 for the Generalized Fermi-Dirac Network (m = 2) and for the Generalized Bose-Einstein Network 
(m = oo) with p = 0.9. The networks have links with energies following a Poisson distribution g{oj) with average c = 10. The 
data are reported for networks of size N = 10,000 (averaged 30 times), N = 5000 (averaged 60 times) N = 2500 (averaged 90 
times) and N = 1250 (averaged 120 times). 
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FIG. 20: (Color online) The maximal distance D from the initial triangle, the maximal degree kmax, and the entropy rate 
Tfll] _|_ h 12 ] 

are plotted as a function of time t for the Generalized Fermi-Dirac Network {m = 2) and for the Generalized 
Bose-Einstein Network (m = oc) with p = 0.9. The inverse temperatures are /3 = 0.05 and /3 = 5 respectively, below and above 
the phase transitions. The networks have nodes with energies following a Poisson distribution g{uj) with average c = 10. The 
data are averaged 20 times. 


In Figure we show major geometrical and structural 
properties of the network as a function of time for given 
values of P below and above the transition for p = 0.9. 
Finally in Figures and we show the degree distri¬ 
bution P{k), the average clustering coefficient C{k) of 
nodes of degree k, and the distribution of the curvature 
P(i?), for the Generalized Fermi-Dirac Network and for 
the Generalized Bose-Einstein Network below and above 
the phase transition. 

In the Generalized Fermi-Dirac Network, for P < f3c^ 
D grows logarithmically with time,( i.e. the network is 
small-world) and the degree distribution P{k) is expo¬ 
nential. For 13 > D grows as a power-law with time, 
(i.e. the network is not any more small-world) k^ax 
increases significantly and P{k) follows a power-law. 


For every value of [3^ the network has high modularity 
M and a hierarchical structure m with an average 
clustering coefficient C{k) of nodes of degree k decaying 
as C{k) k~^ and a = 1. The curvature distribution 
P{R) has a negative tail. 

In the Generalized Bose-Einstein Network, for f3 < f3c^ 
kmax scales sub-linearly with the network size, D in¬ 
creases logarithmically with the network size, (i.e. the 
network is small-world) and the degree distribution P{k) 
is scale-free, the network has high modularity M and 
P{R) has a scale-free positive tail. For /3 > /3c, ^max 
scales linearly with the network size, (i.e. the largest 
node is linked to a finite fraction of the links), D de¬ 
creases significantly, P{k) is dominated by outlier hubs, 
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FIG. 21: Structural and geometrical properties of the Gener¬ 
alized Fermi-Dirac Network as a function of the inverse tem¬ 
perature [3 for single network realizations of size N = 10^. 
The degree distribution P{k), the average clustering coeffi¬ 
cient C{k) of nodes of degree k and the distribution of the 
curvatures P{R) are plotted for /3 = 0.05,0.5,5 and p = 0.9. 
The networks have nodes with energies following a Poisson 
distribution g{uj) with average c = 10. For low values of /3, 
i.e. /3 < /3c, both P{k) and P{R) are exponential, while for 
large values of /3, i.e. > /3c, they become power-law. The 
average clustering coefficient C{k) of nodes of degree k C{k) 
always goes like C{k) oc k~^. 


FIG. 22: Structural and geometrical properties of the Gen¬ 
eralized Bose-Einstein Network as a function of the inverse 
temperature /3 for single network realizations of size N = 10^. 
The degree distribution P{k), the average clustering coeffi¬ 
cient C{k) of nodes of degree k and the distribution of the 
curvatures P(R) are plotted for /3 = 0.05,0.5,1 and p = 0.9. 
The networks have nodes with energies following a Poisson 
distribution g{u}) with average c = 10. For low values of ^d, 
i.e. yd < both P{k) and P{R) are scale-free while for large 
values of /3, i.e. /3 > /3c, these distributions become dominated 
by outliers. The average clustering coefficient C(k) of nodes 
of degree k always goes like C{k) oc k~^ with a <1. 


the network has low modularity M and P{R) has a posi¬ 
tive tail dominated by outliers. For every value of /3, the 
network has a hierarchical structure [78] with an average 
clustering coefficient C{k) of nodes of degree k decaying 
as C{k) k~'^ and a < 1. 


IX. GENERALIZED QUANTUM NETWORK 
EVOLUTION 


A. The evolution 


following, 

IV’Jv(i)) = Ut lipNit - 1)) 

yi — P + ^Jp ^ g-/5(Ci'j'+ej'r')/2 



Wt+2 i,j\i<j 


- 1 ) 
Z(t) 


Here we consider a generalized quantum network state 
evolution corresponding to the generalized geometric net¬ 
work evolution using a similar approach as in the case 
p = 0. We start from an initial condition given by 


IV’Jv(l)) = 


1 




E 




i=l,3 






\i’N{t - 1 )) 


(119) 


where Z{t) is fixed by the normalization condition 
('07v(t)|'0Ar(^)) = 1- Moreover the quantum operator h\-^ 
in Eq. 
state, 


119|) depends on the type of the quantum network 


4243^13 


1001,18) 


where Z{1) enforces the normalization condition 
('0Ar(l)|'0Ar(l)) The evolution of the quantum network 
state is a non-equilibrium Markovian dynamics obtained 
applying the unitary operator Ut to the state |'0Ar(t — 1)). 
This dynamics and the operator Ut are defined as in the 


hL = 


4 


dl 


Generalized Fermi-Dirac 
Quantum Network State 

Generalized Bose-Einstein 
Quantum Network State 


With a long but straightforward calculation, following 
the same steps as for the original Fermi-Dirac Network 








































27 


state and the original Bose-Einstein network state, it is 
possible to show that the normalization constant Z{t) 
can be interpreted as a path integral over the paths cor¬ 
responding to the evolution of the Generalized Growing 
Geometric Networks allowing us to obtain the General¬ 
ized Fermi-Dirac Network for m = 2 and the Generalized 
Bose-Einstein Network for m = oo. 


X. DUAL NETWORKS AND CONNECTION 
BETWEEN THE FERMI-DIRAC QUANTUM 
NETWORK AND SPIN NETWORKS 


Starting from the Fermi-Dirac Quantum network we 
can construct the dual network by mapping each triangle 
of the original network to a node of the dual network and 
every link of the original network to a link in the dual 
network (see Figure 23). 

Each node of the dual of the Fermi-Dirac Network has 
degree 3. A link of the dual network can be saturated or 
unsaturated. A link of the dual network is saturated if it 
connects two nodes of the dual network. This happens if 
and only if the two corresponding triangles of the original 
network are glued together (i.e. they have a common 
link). A link of the dual network is unsaturated if it does 
not connect two nodes. This happens if a triangle of the 
original network has an unsaturated link. 

As the Fermi-Dirac network grows, also the dual net¬ 
work grows. In the original Fermi-Dirac Network only 
process (a) takes place, therefore the dual network is 
a tree. If instead one considers the Generalized Fermi- 
Dirac Network, in which also process (b) takes place, the 
dual network contains loops (see Figure [2^. 

In the case in which the energies of the nodes take 
only integer values we can interpret the dual network as 
a spin network [MSI- In fact we can associate to the link 
^ = {hj) of the dual network the half-integer spin Jij 
given by Eq. <§ satisfying the Clebsch-Gordon condi¬ 
tions given by Eq. ^ at each node of the dual network. 
Nevertheless, we note here that the quantum evolution of 
the proposed Fermi-Dirac network is a non-equilibrium 
dynamics and not an equilibrium one as usually assumed 
in the context of spin networks. 

In the case of the Bose-Einstein Networks it is also pos¬ 
sible to construct the dual network following a similar 
procedure used for constructing the dual of the Fermi- 
Dirac Network, but in this case the dual will not be 
regular, since each triangle in the Bose-Einstein network 
model can be linked to an arbitrarily large number of 
other triangles incident to the links at its boundary. 

It is to mention that in the literature of quantum grav¬ 
ity spin networks with a causal structure have been pro¬ 
posed and are the so called energetic causal sets ga- 
[48] . It will be therefore interesting to explore further 
the connections between the Fermi-Dirac Networks and 
energetic causal sets. 

Our framework is instead quite far from the spin net¬ 
works used in Loop Quantum Gravity. Notably in com- 


Fermi-Dirac 

Geometric 

Network 


Process Processes 

(a) (a) and (b) 
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FIG. 23: (Color online) Fermi-Dirac Network and its dual 
network. In the case in which only process (a) takes place 
the dual network is a tree. In the case in which also process 
(b) takes place the dual network contains loops. In this case 
we have explicitly drawn the triangles forming the underlying 
simplicial complex of the Fermi-Dirac networks. The nodes of 
the dual network are plotted as red squares and correspond to 
the triangles of the Fermi-Dirac network. The links i of the 
dual network correspond to the links {i,j) between nodes i 
and j of the Fermi-Dirac network and are associated to a spin 
variable — Jij = ezj/2. The links of the dual are indicated 
with red lines and can be saturated (if they join two nodes 
of the dual network) or unsaturated if they connect only to a 
single node of the dual network. 


plex quantum network geometries we do not make use 
of intertwines, and the network does not have relevant 
simple symmetries. 


XI. RELATION TO TRIANGULATIONS, 
FOAMS, AND PLANAR GRAPHS 

Planar complex networks have already been studied in 
the literature in several contexts [26ll29l l87ll89] , includ¬ 
ing the study of glass and foams, and planar complex 
networks. 

The model most closely related to our model is the 
scale-free random network constructed by adding ran¬ 
domly triangles to links |26]. Nevertheless, our model 
does not reduce to this model for any value of the param¬ 
eters. In fact also the Bose-Einstein Network for /3 = 0 is 
not equivalent to this model. The difference is that in the 
Bose-Einstein network at /3 = 0, each link is not chosen 
randomly, but proportionally to the number of triangles 
already incident to it (i.e. 1 + according to a kind 
of ’’preferential attachment” to the link. 

Other planar scale-free network models are the pseudo¬ 
fractal scale-free network and the Apollonian net¬ 
works [28]. These network models are deterministic and 
yield scale-free networks with given power-law exponent, 
while the Bose-Einstein Networks are stochastic planar 
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scale-free networks whose power-law exponent depends 
on the distribution g{uj) and on the value of the inverse 
temperature /3. 

Other models for complex networks embedded into 
surfaces have been recently proposed [29] extending ap¬ 
proaches used already in the study of glasses and foams 
[88l l89] . In [29] maximal embedded graphs have been 
characterized using a Monte Carlo algorithm determined 
by an Hamiltonian which is a function of the degree of the 
nodes. This dynamics can explore networks embedded in 
surfaces of different genus, and displays a dynamical slow¬ 
ing down as a function of the inverse temperature of the 
Monte Carlo algorithm. Therefore in these simulations 
the ground state, ordered network is not observed, sim¬ 
ilarly to what has been observed in models of glass and 
foam dynamics [88] |89] . This approach is very different 
from the one proposed in the present article, although 
the focus is always the characterization of the network 
geometry. For example the phase transitions present in 
the Quantum Complex Network Geometries are very dif¬ 
ferent from the one observed in [29] |88j [89] . In fact in 
the Complex Quantum Network Geometries the observed 
phase transitions are non-equilibrium phase transitions, 
they are determined by the quenched disorder in the net¬ 
work. Instead, in [29] HSl [89] the network has an equi¬ 
librium Hamiltonian dynamics, the ground state is well 
defined but there is a dynamical slow down. Moreover, 
in Gomplex Quantum Network Geometries the observed 
phase transitions can change the metric properties of the 
networks, which can go from a small-world network to a 
network with large diameter without changing its Euler 
characteristic, as in the case of the phase transition ob¬ 
served for Fermi-Dirac networks. Instead in [29], changes 
from small world networks to networks with large diam¬ 
eter occur as a function of the Euler characteristic of the 
network. 


XII. CONCLUSIONS 

In this work we have proposed to study a geometri¬ 
cal network evolution based on growing simplicial com¬ 
plexes of dimension 2, i.e. simplicial complexes formed 
by triangles. The network evolution describes the evo¬ 
lution of a quantum network state defined by a path in¬ 
tegral. The quantum network states are characterized 
by a set of quantum occupation numbers that can be 
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